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1. Introduction 

Nonlocal spatial solitons have been investigated for decades both theoretically and exper- 
imentally. Spatially nonlocality is a generic property in different materials including the 
photorefractive media [1, 2], thermal nonlinear media [3-7], liquid crystals [8-13], and so 
on. Nonlocality can lead to new kinds of waves that would have been otherwise impossible 
in local nonlinear media. Especially, for two-dimensional media with different types of non- 
local response, they can support stationary multipoles [14-19], stable vortices [20-24], and 
rotating [25, 26] and spirahng [5, 27] soliton states. 

Surface waves localized at the interface of two different optical materials have many novel 
properties, which have been studied in nonlocal media recently. The light beam trajectory 
can be strongly affected by the presence of interfaces, because beams propagating in nonlocal 
media cause refractive index changes in regions far exceeding the beam width. Under proper 
conditions, stationary surface waves can propagate along the interface in both local nonlinear 
media [28-31] and nonlocal nonhnear media [32-39]. 

Thermal media [34, 35, 38] and photorefractive crystals [39] have been also utilized to 
demonstrate (2+l)D surface solitons. In this paper, we will study (2+l)D surface solitons 
occurring at the interface between a linear medium and a nonlocal nonlinear medium whose 
nonlinear contribution to the refractive index has a initial value at the interface. We find 
that there exist single surface solitons at the edge and the corner of the two 2D media and 
dipole surface solitons at the edge of the two 2D media. These stable solitons do not exhibit a 
power threshold. The positions of the peak values and full width at half maximum(FWHM) 
of the surface solitons can be affected by the degree of nonlocality. However, the initial value 
can only influence the positions of the peak values of the surface solitons. In addition, when 
a laser beam is launched away from the interface, the beam will be periodic oscillations, 
even if the launch position is far away from the interface. 

2. Theoretical Model 

Considering a laser beam propagating along the interface between a nonlocal nonlinear 
medium and a linear medium, the complex amplitude E{X^ Y, Z) of the light field satisfies 
the scalar wave equation [32, 34, 39] 

^'xyE + ^ + kyE^O, (1) 
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where fco = 27r/A is the wave number in vacuum, and n — ul + An for the nonhnear 
medium(X < 0) and n — no for the hnear medium(X > 0). An represents the nonhnear 
contribution to the refractive index and may originate from any diffusive nonhnear effect, 
which can be written by [32] 

t(;^V^yAn - An + n2\E\^ = 0, (2) 

where Wm is the characteristic length of the nonhnear response and n2 is the nonhnear index 
coefficient. For the local case, Wm 0, we have An = n2\E\'^. 

Let us put E{X^Y^Z) = £'o(X, F, Z) exp(i/3Z) and submit into the equation (1) and 
(2). Then, using the slowly varying envelope approximation and introducing the normalized 
variables x — X/wq^ y — Y/wq^ z — Z/ {/3wq)^ a — k()W()^Jn^nLEQ and (j) — k^w^nLAn^ we get 

^ V^a + idza + /3ia + 0a = for x < 0, (3a) 

-Via + id^a + /^sa = for a; > 0, (3b) 
2 

and 

rf^Vi0-(/)+|ap = O for x<0, (4) 

where = 5^ + 5^, /?i = w\{kl^n\^ - /?^) /2, /?2 = w\{k\r?o - 0^) /2, /? is the wave number in 
the media, w^) is the beam width, and d — Wm/wQ stands for the degree of nonlocality of the 
nonlinear response. For x < 0, the equations describe a local nonlinear response as g? ^ 
and a strongly nonlocal response as d ^ oc. 

We search for stationary soliton solutions of Eqs. (3) and (4) numerically in the form 
a(x, z) — u{x^ y) exp{ibz)^ where u is the real function and 6 is a real propagation constant 
of spatial solitons in the normalized system. 

-V^^u -bu + f3iu + (j)u^O for x < 0, (5a) 
2 

^ V^M -bu + ^2U = for a; > 0, (5b) 

and 

ImI^ = for a;<0. (6) 
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3. Numerical results 

A. Single surface soUtons 

We firstly talk about new surface-wave soliton solutions at the edge of the two 2D media. 
Here, we assume that the normalized nonlinear contribution to the refractive index at the 
interface(y = 0) has a initial value 0^(0^ > 0). The boundary conditions for the fields at 
the interface are the continuity of the transverse field(a(y ^ 0_) = a{y 0+)) and its 
derivative(G?a(?/ 0-)/dy = da{y 0^)/dy). Because the width of the surface solitons is 
much smaller than the sample width, a and (/) vanishes at the other boundary. 

Three diflFerent solitons are separately shown in Fig. 1(a), (b) and (c). We can easily see 
that the soliton is closely attached to the interface with (l)d = 7 and d = 20 in Fig. 1(b), but 
the solitons are farther detached from the interface with (j)^ — 1 and = 20 in Fig. 1(a) or 
(t)d — 7 and = 5 in Fig. 1(c). Obviously, and d influence the position or shape of the 
solitons. To further explain this point, we see that the positions of the peak values (?/^a^c) 
and FWHM of the surface wave solitons versus the boundary value (/)^(Fig. 1(d)) at the 
interface or the degree of nonlocality rf(Fig. 1(e)). From Fig. 1(d), we can see that, under 
the condition((i = 20), the soliton will be attracted to the interface and more and more 
significant part of their optical power residing in the linear meidium as (j)^ increases. That 
is to say, is larger, the larger a "surface force" exerted on the beam by the interface. 
However, 0^ cannot influence the beam width of solitons. From Fig. 1(e), at = 7, the 
changing of the refractive index of the nonlocal nonlinear medium induced by d results in 
the changing of FWHM and ymax of the solitons. Because of the changing of FWHM, ymax 
changes intricately, though a force exerted on the beam by the degree of nonlocality increases 
all the while. 

Fig. 2(a) shows that the energy flow U — /^/^|aprfxrfy of the single 2D surface 
solitons monotonically increases with h where dU/db > 0. This shows that the solitons are 
stable [33, 34]. Here, to further elucidate the stability of the surface solitons, we do the 
direct numerical simulations of Eqs. (3) and (4) with input conditions a\z=o = u{l + p), 
where p{x^y) is a broadband random perturbation. The fact which is shown in Fig. 2(b) 
confirms the result of Fig. 2(a). Then we proceed to address the dynamics behavior of the 
propagation of surface solitons. For convenience, the (1+1)D circumstance is considered. 
Fig. 2(c) depicts that a narrow beam is launched y — — 1.99/im away from the interface. The 
beam maintains a localized shape. However, it is oscillation in a fully periodic fashion in the 
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virtue of the cooperation of the forces exerted by the boundary and the nonlocal nonlinearity, 
even if the launch position is far away from the opsition of the surface soliton [34] . 

Next, we consider that new soliton solutions at the corner of the two 2D media. Here, 
we assume that the normalized nonlinear contribution to the refractive index at x = and 
y = have a initial value (j) = (l)d{(l)d > 0). The boundary conditions for the fields at the 
interface meet the continuity conditions. 

Fig. 3(a), (b) and (c) separately show three different solitons. The soliton is closely 
attached to the interface with (j)^ = 10 and d = 20(see Fig. 3(b)), but the solitons are 
farther detached from the interface with 0^ = 5 and d — 20(Fig. 3(a)) or (/)^ = 10 and d — 
10(Fig. 3(c)). To further illustrate the influence of (j)d and d on the solitons, we display that 
ymax and FWHM versus (j)d in Fig. 3(d) or d in Fig. 3(e). The energy flow U monotonically 
increasing with b shown in Fig. 4(a) and the direct numerical simulations of Eqs. (3) and (4) 
with noise cr'^oise — ^-^^ shown in Fig. 4(b) explain that the solitons are stable. The results 
can be similarly illustrated as the solitons at the edge of the interface. 
B. Two dimensions dipole surface solitons 

In addition to single surface solitons, we also flnd a 2D stationary dipole surface solitons. 
The surface solitons are found numerically by a standard relaxation method which converges 
to a stationary solution after some iterations provided that a suitable guess for initial field 
distribution. The boundary conditions are the same as the single surface solitons at the 
edge of the interface. 

Fig. 5(b) depicts that the amplitude for a dipole soliton which is closely attached to the 
interface with (j)^ = 10 and d = 20. However, the solitons are farther detached from the 
interface with (j)d — ^ and d = 20 in Fig. 5(a) or 0^ = 10 and d — 10 in Fig. 5(c). Because the 
poles of solitons are almost symmetric in the x direction, the changing of ymax and FWHM 
of one pole of the solitons can explain the changing of the position and shape of the dipoles. 
Fig. 5(d) and (e) shows ymax and FWHM as functions of the boundary value and the 
degree of nonlocality respectively. 

Fig. 6(a) shows that the energy flow monotonically increases with the propagation con- 
stant for dipole surface solitons. With increasing energy flow, surface dipoles become more 
localized, i.e., the distance between poles along the x axis and their widths decrease. To 
further elucidate the stability of the surface dipoles, we do the direct numerical simulations 
of Eqs. (3) and (4) with noise cr'^oise ^ O.Ob. Particularly, the complex surface solitons are 
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stable in the entire existence domain. Fig. 6(b) depicts a typical evolution dynamic. The 
considerable input perturbations cannot almost cause oscillations of amplitudes of the two 
poles forming the dipole, but the dipoles remain their internal structures over huge distances. 
Contrarily, in a bulk diffusive medium, it is known that weak input perturbations can cause 
slow but progressively increasing oscillations of the bright spots forming a dipole, resulting 
specially in their slow decay into fundamental solitons [35]. So, we think that the presence 
of a interface possessing a initial value for </) leads to stabilization of dipole solitons. This is 
further illustrated by the results of Fig. 6(c). when the beam is launched y = — 1.91/im away 
from the interface, the stationary soliton cannot be formed, but it is periodically oscillating. 
This are very good description of the boundaries of the role of solitons. 

4. Conclusion 

In conclusion, we have addressed (2+l)D surface solitons occurring at the interface between 
a linear medium and a nonlocal nonlinear medium whose nonlinear contribution to the 
refractive index has a initial value at the interface. We find that there exist stable single 
and dipole surface solitons which do not exhibit a power threshold. The degree of nonlocality 
have influence on the positions of the peak value and FWHM of the surface solitons, but the 
initial value can only influence the positions of the peak values of the surface solitons. In 
addition, when a laser beam is launched away from the interface, the beam will be periodic 
oscillations, even if the launch position is far away from the interface. 
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List of Figure Captions 

Fig. 1. Sketch of 2D single surface solitons at the edge of the interface with (a) 0^ = 1, 
d = 20, (b) 0^ = 7, = 20 and (c) 0^ = 7, = 5. The positions of the peak values i/max and 
FWHM versus (d) the boundary value 0^ and (e) the nonlocal degree d. White dashed line 
indicates interface position. All quantities are plotted in arbitrary dimensionless units. 
Fig. 2. (a) Energy flow U versus the propagation constant b with g? = 20 and (j)d = 7. (b) 
Stable propagation of surface solitons in Fig. 1(b) with noise cr'^oise — O-O^ ^ distance of 15 
diffraction lengths, (c) Trajectories of the incident beam with the beam center coordinates 
y = — 1.99/im. White dashed line indicates interface position. All quantities are plotted in 
arbitrary dimensionless units. 

Fig. 3. Sketch of 2D single surface solitons at the corner of the interface with (a) = 5, 
= 20, (b) = 10, = 20 and (c) = 10, = 10. The positions of the peak values ymax 
and FWHM versus (d) the boundary value (/)d and (e) the nonlocal degree d. White dashed 
line indicates interface position. All quantities are plotted in arbitrary dimensionless units. 
Fig. 4. (a) Energy flow U versus the propagation constant b with = 20 and (/)^ = 10. 
(b) Stable propagation of surface solitons in Fig. 3(b) with noise (Jnoise ^ O-O^ for a distance 
of 15 diffraction lengths. White dashed line indicates interface position. All quantities are 
plotted in arbitrary dimensionless units. 

Fig. 5. Sketch of 2D dipole surface solitons with (a) (l)d — 1, = 20, (b) (pd — 10, g? = 20 
and (c) = 10, G? = 10. The positions of the peak values ymax and FWHM versus (d) 
the boundary value and (e) the nonlocal degree d. White dashed line indicates interface 
position. All quantities are plotted in arbitrary dimensionless units. 

Fig. 6. (a) Energy flow U versus the propagation constant b with d = 20 and (/)d — 10. (b) 
Stable propagation of surface solitons in Fig.5(b)with noise cr'^oise — O-O^ ^ distance of 15 
diffraction lengths, (c) Trajectories of the incident beam with the beam center coordinates 
y — — 1.91/im. White dashed line indicates interface position. All quantities are plotted in 
arbitrary dimensionless units. 
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Fig. 1. Sketch of 2D single surface solitons at the edge of the interface with (a) 0(^ = 1, 
d = 20, (b) (/)d = 7 ^ d = 20 and (c) (pd = 7^ d = 5. The positions of the peak values i/max and 
FWHM versus (d) the boundary value (/)d and (e) the nonlocal degree d. White dashed line 
indicates interface position. All quantities are plotted in arbitrary dimensionless units. 
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Fig. 2. (a) Energy flow U versus the propagation constant b with d = 20 and (pd = 7. (b) 
Stable propagation of surface sohtons in Fig. 1(b) with noise (Jnoise = O-O^ for a distance of 15 
diffraction lengths, (c) Trajectories of the incident beam with the beam center coordinates 
y = — 1.99/xm. White dashed line indicates interface position. All quantities are plotted in 
arbitrary dimensionless units. 
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Fig. 3. Sketch of 2D single surface solitons at the corner of the interface with (a) (j)d = 5, 
d = 20, (b) (/)d = 10, = 20 and (c) (pd = 10, = 10. The positions of the peak values i/max 
and FWHM versus (d) the boundary value (/)d and (e) the nonlocal degree d. White dashed 
line indicates interface position. All quantities are plotted in arbitrary dimensionless units. 
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(b) 



Fig. 4. (a) Energy flow U versus the propagation constant b with d = 20 and (pd = 10. 
(b) Stable propagation of surface sohtons in Fig. 3(b) with noise (Jnoise = O-O^ for a distance 
of 15 diffraction lengths. White dashed line indicates interface position. All quantities are 
plotted in arbitrary dimensionless units. 
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Fig. 5. Sketch of 2D dipole surface solitons with (a) = 1, = 20, (b) (j)d = 10, = 20 
and (c) = 10, = 10. The positions of the peak values ymax and FWHM versus (d) 
the boundary value (^a and (e) the nonlocal degree d. White dashed line indicates interface 
position. All quantities are plotted in arbitrary dimensionless units. 
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Fig. 6. (a) Energy flow U versus the propagation constant b with d = 20 and (pd = 10. (b) 
Stable propagation of surface sohtons in Fig.5(b)with noise (Jnoise = O-O^ for a distance of 15 
diffraction lengths, (c) Trajectories of the incident beam with the beam center coordinates 
y = —1.91fim. White dashed line indicates interface position. All quantities are plotted in 
arbitrary dimensionless units. 
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